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ABSTRACT 



O The Cabibbo suppressed decays D — » V7 + Z (V is light vector meson) 

present in principle the opportunity to observe the short distance FCNC tran- 
sition c — > ul + l~ , which is sensitive to physics beyond the Standard Model. 
We analyze these as well as the Cabibbo allowed D — > Vl + l~ decays within 
the Standard Model, where in addition to the short distance dynamics also the 
long distance dynamics is present. The long distance contribution is induced 
by the effective nonleptonic weak Lagrangian accompanied by the emission of a 
virtual photon, which occurs resonantly via conversion from a vector meson p°, 
lo or (ft or nonresonantly as direct emission from a D meson. We calculate the 
branching ratios for all D — > Vl + l~ decays using the model, which combines 
heavy quark symmetry and chiral perturbation theory. The short distance con- 
tribution due to c — > ul + l~ transition, which is present only in the Cabibbo 
suppressed decays, is found to be three orders of magnitude smaller than the 
long distance contribution. The branching ratios well above 10 -7 for Cabibbo 
suppressed decays could signal new physics. The most frequent decays are the 
Cabibbo allowed decays, which are expected at the rates, that are not much 
lower than the present experimental upper limit: Df — > p + /j, + fj>~ is expected at 
the branching ratio of approximately 3- 10 -5 , while D° — » K*°[i + ii~ is expected 
at 1.7- 1(T 6 . 



I. INTRODUCTION 



In the charm sector phenomena like D° — D° mixing, CP-violation and rare 
decay probabilities are small, which makes them good candidates as probes 
for new physics with small background from the Standard Model ||, ||, H]. In 
particular, decays of type D — ► Xl + l~ were sing led out g |, | as a possi- 
ble good window to non-standard contributions of the flavour-changing neutral 
transition (FCNC) c —* ul + l~~ , at the 10 level for the branching ratios. This 
suggestion was prompted by the smallness of the short-distance (SD) c — > ul + l~ 
contribution within the Standard Model, which leads j|, ||] to a branching ra- 
tio of only 10 -9 for the inclusive process. Although QCD corrections to this 
process have not been calculated in detail yet, these are not expected to affect 
significantly the size of the c — ► ul + l~ amplitude, as explained in the next sec- 
tion. Accordingly, one expects the hadronic exclusive decays induced by this 
SD transition to occur with branching ratios of the order of 10~ 10 . 

Further studies, which have considered the long-distance (LD) contribu- 
tion to D — * Pl + l~ transitions (P is light pseudoscalar) || fjj have concluded 
that these are larger than SD ones. The analysis of the LD contributions in 
D + >° — > ir + '°l + l~ @ has shown these modes are expected to lead to branching 
ratios of the order of 10 -6 in the resonance region and a few times 10~ 7 in 
the nonresonant region, thus practically invalidating their use for observing the 
c — > ul + l~ transition within the Standard model. 

A similar situation holds in the case of D — > Vj (D — * P7 is forbidden), 
where long distance effects [§, [f| 10] cause these modes to have branching ra- 
tios in the 10~ 7 — 10~ 4 range. On the other hand, the SD component due to 
the magnetic electroweak penguin transition c — > wy is GIM suppressed and 
does not reach beyond the 10~ 9 range at most, despite of being considerably 
enhanced by QCD corrections H, O]. Thus, here again the LD effects mask 
the contribution of the SD c — > wy loop, except for very special circumstances, 
which were pointed out recently || 12 1. 

We analyze LD and SD contributions to all D — > Vl + l~ decays within 
the Standard Model. The SD contribution due to c — * ul + l~ is present only 
in the Cabibbo suppressed decays D° — ► p°l + l~, D° — ► luI + 1~ D° — > cj)l + l~, 
D + — > p + l + l~ and Df — > K* + l + l~. Our results should provide the appro- 
priate theoretical background against which possible signals of new physics 
are searched for in these decays. Motivated by the experimental searches, we 
analyze also the Cabibbo allowed decays (D° — ► K*°l + l~ and Df — ► p + l + l~), 
which are the best candidates for their early detection, and the doubly Cabibbo 
suppressed decays (D + — ► K* + l + l~ and D° — ► K*°l + l~). Here the signals from 
new physics are not expected from the theoretical models usually considered. 
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On the experimental side, so far there are only upper bounds on the branch- 
ing ratios of D — > Vl + l~ decays from E653 and CLEO [13, 14], in the range 
1CT 3 — 10~ 4 , but these are expected to improve in the future. 

In Sec. II we present the details of our approach and we define the ap- 
proximations used. In Sec. Ill we give the results of our calculations and we 
summarize in Sec. IV. 



II. MODEL DESCRIPTION 



A. Long distance contributions 



In this subsection we present the general framework used for calculating the 
long distance amplitudes, while the details of the model employed are given in 
subsection II C. The long distance contribution in D — > Vl + l~ decays is due to 
the effective nonleptonic weak Lagrangian, which induces the weak transition 
between the initial and final hadronic state. The weak transition has to be 
accompanied by the emission of a virtual photon, which finally decays into a 
lepton antilepton pair. The effective nonleptonic weak Lagrangian responsible 
for charm meson decays is 

G F 

£ld = --y=V uqi V^. [a^uqiYiqjc)^ + a 2 (uc) fl (q j q i ) fl } : (1) 

where (V'iV'2) M = V'lT^U — T 5 )^, Qi,j represent the fields of d or s quarks, Vij 
are the CKM matrix elements and Gf is the Fermi constant. In our calculation 
we use a\ = 1.26 and a 2 = —0.55 as found in Jl5| , from an extensive application 
of (H) to the study of nonleptonic D decays. 

The virtual photon emission from the hadronic states is taken in our ap- 
proach to proceed through two different mechanisms: 

(i) In the nonresonant mechanism the photon is emitted directly from the initial 
D state. 

(ii) In the resonant mechanism, apart from the final vector meson V, an ad- 
ditional neutral vector meson Vq is produced, which converts to a photon 
through vector meson dominance (VMD). In this case, a nonleptonic weak 
decay D — > Wq is followed by the transition Vq — ► 7* — > l + l~, where Vq is a 
short-lived vector meson p°, uj or (p. 

The evaluation of the matrix elements of the product of two currents (||) 
requires nonperturbative techniques and we are forced to use some approxima- 
tion. We have undertaken to use systematically the factorization approxima- 
tion, where the matrix element of the product of two currents is approximated 
by 

{v^MiYiqkc)^) = (vmyrm^q-kc)^) 
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+ (l\(q iqj Y\V)(V\(q k c),\D) 

+ (Vj\(q iqj r\0)(0\(q k c)„\D) . (2) 

The first two terms are the spectator contributions, in the following denoted by 
As P ec,-y and As pe c,v, respectively, and the third term is the weak annihilation 
contribution, denoted by AAnnih- Here 7 denotes the virtual photon. 

To calculate the matrix elements in (pi) we use the hybrid model, which 



combines heavy quark effective theory (HQET) and chiral Lagrangian |16||-[18|, 
which has been successfully employed already for D meson decays in several 
papers |l^]-[^6|. A very detailed description of the hybrid model and its previous 
applications is given in fl8|| . 

The relevant hadronic degrees of freedom for the present calculation are 
heavy pseudoscalar (D) and vector mesons (D*) and light pseudoscalar (P) 
and vector (V) mesons. Within this approach the diagrams that contribute 
to the amplitudes As pec> j, As pec ,v and AAnnih (@) are shown in Figs, la, lb 
and lc, respectively. Different diagrams in Fig. 1 are denoted by the roman 
numbers from / — VIII. The diagrams III and IV represent nonresonant con- 
tribution (mechanism (i)). All the remaining diagrams, which proceed through 
the intermediate short-lived vector meson Vq (p, to and </>), represent the reso- 
nant contribution (mechanism (h)). The resonant amplitude is represented in 
the whole q 2 region by the Breit-Wigner vector meson propagator, (q is the 
sum of lepton and antilepton momenta). In the regions of q 2 far away from 
rnyj, the resonant amplitude is given therefore solely by the tail of the Breit- 
Wigner vector meson propagator. The square in each diagram of Fig. 1 denotes 
the weak transition due to the effective Lagrangian Cld ([]])• This Lagrangian 
contains a product of two left handed quark currents (qkqi) 11 , each denoted by 
a dot on Fig. 1. The left handed currents will be expressed in terms of the 
relevant hadronic degrees of freedom: D, D* , P and V. In our notation the 
hadronic current J2 in the diagram II, for example, creates V meson, while the 
hadronic current J\ annihilates D and creates Vo at the same time. 

In the model we use, there is no contribution of J/\I/ or other cc excited 
states. The contribution that would arise by the exchange of this mesons is ef- 
fectively described by the diagram III of Fig. 1, where cc exchange is "hidden" 
in the DD*^ coupling. The alternative approach of a direct cc exchange would 
require the knowledge of their couplings to photons over a wide region of q 2 , of 



which one has only rudimentary knowledge [27]. 



B. Short distance contributions due to c — ► ul + l 







In addition to long distance dynamics, the Cabibbo suppressed decays D 
p°l + l~, D° -> wl+r, D° -► </>l+l-, D+ -» p+l+r and D+ -» K*+l+l~ can 
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also be driven by the short distance c — ► ul + l~ transition. The short distance 
part in D — > Vl + l~ decays will turn out in general to be much smaller than the 
long distance part. However, we shall find that in the case of D° — ► p°(uj)l + l~ 
the short distance part is of the same order of magnitude as the nonresonant 
part of the long distance contribution. In this subsection we estimate the size 
of the short distance amplitudes, which is nonresonant in its nature. 

The effective Lagrangian for FCNC transition c — ► ul + l~ arises from WW 
exchange box diagrams and Z and 7* penguin operators 0. It has been ob- 



tained using the similar results for s — > dl + l decay [28 
G F e 2 



C 



SD 



y/2 I67r 2 sin 2 e w 

Vi n 7At (l- 75 )c UiT-fO— ys)l + BiT-f(l 

i=d,s,b 



(3) 



where the Willson coefficients Ai, Bi and F\ are given in Appendix A and Vi 
are the CKM coefficients, Vi = V£V u i. The expression @ does not contain 
the QCD corrections, which have not been studied for c — * ul + l~ decays so 
far. When referring to these corrections for c — > ul + l~ , one is reminded that 



in the case of c — > wy decay there is a huge QCD enhancement |8|, 11], which 
is due to the following reason: The effect of QCD is that the Wilson coeffi- 
cient cj{m c ), responsible for the magnetic penguin decay c — > wy, obtains the 
admixture of the other Wilson coefficients evaluated at the scale m\y in the 
leading order Ci(myy), i = 1..10. Since cy(m^) is extremely suppressed com- 
pared to some other Wilson coefficients Ci(m\v), the resulting cj(m c ) is much 
bigger than ci(m\v) and a huge QCD enhancement for c — > wy occurs. In the 
c — ► ul + l~ decay on the other hand, the responsible Wilson coefficients A(mw) 
and B(myy), which do not contribute for the real photons, are not suppressed 
in the lowest order and one would not expect large QCD effects on them as one 



has learned from estimation of s — > dl + l~ [29]. The coefficient of the magnetic 
transition ^(mw), which is proportional to Wilson coefficient cj(mw), indeed 
acquires large QCD correction, however it is strongly suppressed in the lowest 
order. Since we are only interested in the rough estimation of the short distance 
contribution in D — > Vl + l~ and the last term in @ is much less important than 
the other two for such decays p£^| , we neglect the last term altogether, using 
the approximation explained in Appendix A. 

The Lagrangian Csd (§) then gives the branching ratio for inclusive c — > 
ul + l~ process 



T(c^ul+r) G 2 F m\ ( a \ 2 



T(D°) " 1927r 3 r(L>°) V 4tt sin 2 W 



ViAA' + lViBA 2 } =2.9 10" 



-9 
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To predict the exclusive amplitudes for D — > V7 + Z induced by Csd we 
have to evaluate the matrix elements 

(V\u ltl (l-j 5 )c\D) . (4) 

We shall do this by again using the hybrid model, which is described in the 
next subsection. The corresponding Feynman diagrams within this approach 
are given in Fig. 2. The squares in the diagrams denote the weak transition due 
to the short distance Lagrangian Csd ©• This Lagrangian contains a product 
of a quark and lepton weak currents, each denoted by a dot in Fig. 2. We re- 
mark, that these diagrams have a long distance counterpart, given by diagrams 
V and VI of Fig. 1, which represent the long distance c — > wy transition ^ 



C. Theoretical framework: chiral Lagrangians, heavy quark limit 
and vector meson dominance 



Here we present the model, which we use to evaluate the matrix elements 
(||) and (||) needed to predict D — * Vl + l~ amplitudes. The framework we use 
for our treatment is that of an effective Lagrangian, which embodies two impor- 
tant approximate symmetries of QCD, the infinite heavy quark Q mass limit 
(rriQ — ► 00) and the chiral limit for light quarks, namely (m u , m^, m s ) — ► 0. 
This approach, which was developed during the last few years (fll6fl-|26|1 and 
additional references quoted in |Tj|), has been used with a good measure of 
success to treat strong, electromagnetic and weak decays of D and B mesons. 
Obviously, an effective Lagrangian approach has also its weakness, as it involves 
a number of unknown coupling constants. Fortunately, the use of observed pro- 
cesses makes it possible to determine a good proportion of them, as detailed in 
this and the next subsection. Moreover, the use of form factors alleviates the 
limitations on the range in which the basic assumptions of the model hold to 
a good accuracy. The model and its various applications till now are well ex- 
posed in a recent review [|TJ]. In the present subsection, we describe those parts 
which are needed for our calculation. We introduce now the strong and electro- 
magnetic interaction Lagrangians for the heavy (hadrons containing c quark) 
and light (hadrons containing only light u, d and s quarks) sector and the rel- 
evant weak currents. At the end of the section we discuss the values of free 
parameters, that enter the Lagrangians and currents in our model. Our strong 
and electromagnetic Lagrangian |16[|-||18[| is invariant under heavy quark spin 
(SU(2)), chiral (SU(3)l x SU(3)r), Loretz, parity and U(l) gauge transforma- 
tion |TJ] . The light vector mesons are incorporated using the hidden symmetry 
approach [18, [30 1. We are aware that using HQET, which converges very slowly 
in the case of c quark, presents a rather rough approximation. In spite of that, 



5 



the HQET approach, which helps to reduce the number of free parameters, has 
been successfully applied in many D decays (e.g. fTJI and references therein). 
The light degrees of freedom are described by the 3x3 Hermitian matrices 



and 



n 



7T U , T? 8 I rjo 



V2 



7T 



I V8 | VP 
^ V6^ V3 



K+ 
' 7e ~ l ~ Vz' 



(5) 



Pi 

V K*; 



K*° 



K*° 



F^{p) = d^p u - duPf, + [p,j,, p v \ (6) 



/ 



for the pseudoscalar and vector mesons, respectively. They are usually ex- 
pressed through the combinations 



exp(y) 



(7) 



where / ~ f n = 132 MeV is the pion pseudoscalar decay constant and 

P» = l ^7xPn > ( 8 ) 



■ 9V 



where gy is fixed in the case of the exact flavor symmetry to be the VPP 
coupling gy = 5.9 Q. 

The most general strong Lagrangian for the light mesons in the leading 
order of chiral perturbation theory is [BT1] 



r 1 



1 



light ~ ~ J -{tr{A^) + atr[(V, - p^) 2 ]} + ^-tr[F, u (p)F^ (p)] , (9) 



(10) 



where we have introduced two currents 

Vfj, = ^{u^D^u + uDpui) and = -(u^D^u - uD^v)) 

Demanding the Lagrangian @ to be invariant under the local gauge transfor- 
mation, corresponding to the electro-magnetic U(l) transformation in QCD, 
we define the covariant derivatives as 



D^u = (dp + B^u and D^v) = (<9 M + B^i 



t 



with = ieB^Q, Q = diag(2/3, —1/3, —1/3) and B^ being the photon field. 
The constant a @ is in principle a free parameter. We fix it to a = 2 p0|| 
assuming the exact vector meson dominance, where the pseudoscalars interact 
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with the photon only through vector mesons. With this choice, the photon - 
vector meson interaction given by the second term of Lagrangian @ is 

C Vl = -eg v f 2 B,{p^ + ^ - ^V) . (11) 

Instead of using the exact SU(3) symmetry values gy = 5.9 and / = 132 MeV, 
we express the V'j couplings in terms of the measurable quantities g p , g^ and 
g^ defined by the matrix element of the corresponding vector current Jy 

(V(e v , q )\J v \0)=g v (q 2 )e*»( q ) . (12) 

In our calculation we use the values gy given in Table 2, which have been di- 
rectly measured in the leptonic V — > l + l~ decays, and we make the assumption 
9v(q 2 ) = 9v( m y) = 9V- The photon - vector meson interaction Lagrangian 
(|ll|) defined through the constants gy is 

C Vl = -^(9 P P°» + f ^ - ^f>r) Bp. (13) 

As far as the calculation of the amplitudes for the diagrams of Fig. 1 is con- 
cerned, the Lagrangian Cu g ht @ provides also the VVV vertex given by the 
third term in (Q). The VVV vertex is present in the diagram VIII of Fig. 1, 
which describes the photon emission from the charged vector meson. 

We need also the PVV vertex, which is present in the diagram VII of Fig. 
1. This interaction term can be generated only in the next-to- leading order of 
chiral perturbation theory as [18] 



flight = -^j^e^Tr{d,p v d a p^) , (14) 



where Cyyu is free parameter. 

Both the heavy pseudoscalar and the heavy vector mesons are incorporated 
in a 4 x 4 matrix 

Ha = l -^{P* a ^ ~ Pal,) , H a = f = {P$f + ^Ts)^ , (15) 

where a = 1, 2, 3 is the SU(3)y index of the light flavors and P*^, P a annihilate 
a spin 1 and spin heavy meson Qq a of velocity v, respectively. The strong 
and electromagnetic Lagrangian in the heavy sector have to provide us with 
the DDj, DD*j and DD*V vertices. The first vertex describes the photon 
emission from the charged D meson and is generated in the leading order of 
HQET (invariant under heavy quark symmetry and U(l) gauge transformation 
with minimal number of derivatives) as [17, ll^] 



Zhavy = ^r[F a ^(^ + V, - jieBjHj (16) 
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The DD*^ and DD*V vertices can be generated only in the next-to-leading 
order of the heavy quark and chiral expansion and are described by [17, |l^] 



CLavy = -X'Tr[H a a^F^(B)H a ] + i\Tr[H a a^F^ {p) ab H b \ . (17) 

The first term contributes to the diagram /// Fig. 1, while the second term 
contributes to the diagrams / and V of Fig. 1 and the diagram / of Fig. 2. 
The A and A' are free parameters. 

In addition to the strong and electromagnetic interaction, we have to specify 
the weak one. The effective weak Lagrangian responsible for the long distance 
contribution is given by Cld @ and for the short distance contribution by Csd 
(||). As we deal with the probabilities for the weak decays of hadrons, we rewrite 
the quark weak currents in Cld ([[]) an d C-sd (@) in terms of hadronic degrees 
of freedom. The weak current <j a 7^(l — 7 5 )c containing a c quark and one 
light anti-quark q a transforms under chiral SU(3)l X SU(3)r transformation 
as (3l, 1_r)- At the hadronic level we impose the same chiral transformation 
and we require the current to be linear in the heavy meson fields D a and D* a 
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J£ = l -iaTr[^{l-^)H b u\ a ] (18) 
+ ai Tr[ l5 H b (p» - V%utj + a 2 Tr[^ l5 H b v a (p a - V a ) bc u\ a ] + ... . 



The current ([18]) is the most general one in the leading 1 / m c order of HQET and 
next to leading order of chiral perturbation theory. The first term is connected 
to the definition of the heavy meson decay constant (D(p)\(q a c)^\0) = —i/dP^, 
where a = fo^J^D- The second and third term contribute to the diagrams II 
and VI of Fig. 1 and diagram iT of Fig. 2, where J£ (|l8| ) annihilates D meson 
and creates V or Vo meson at the same time. The constants a± and oli are free 
parameters. 

D. The choice of the parameters 

We now turn to the values of the coupling constants Cyvu, A, A', a\ and 
«2, which we need in the evaluation of the amplitudes of diagrams on Figs. 1 
and 2. 

The coupling Cyvu can be determined in the case of the exact SU(3) flavor 



symmetry following the hidden symmetry approach of [30] and is found to be 
[Cyv-nl = 3<^/32vr 2 = 0.33. Experimentally, it can be directly determined from 
the V — > PVq — > Pj decay rates. In the following we will use the average value 
of Cwxii obtained from the measurements of different V — > PVq — > P7 decays 
|CVvn| = 0.31 §§, which is close to its SU(3) limit. 
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We determine the three parameters A, a\ and «2 using three values related to 
the helicity amplitudes Br = 0.048±0.004, T L /T T = 1.23±0.13 and T + /r_ = 
0.16 ± 0.04 for the process 

D + _> K*°l+vi ||], taken from the average of data 



from different experiments [14|. We get four sets of solutions for A, a,\ and 
«2 0] and we choose the set, which gives the best fit with a number of the 
nonleptonic decays D —* PV, D —* VV and D —* PP Q: 
A = -0.34 ± 0.07 , ai = -0.14 ±0.01 , and a 2 = -0.83 ± 0.4 . 

In order to gain information on A' we turn to an analysis of D*° — ► D°j, 
D* + — > D + 7 and D* + — > -D^7 decays. Experimentally, only the ratios i?S^ = 

r(D*° _» D° 7 )/r(L>*° -» L>°vr ) and i?+ = r(z?*+ -> D+ 7 )/r( J D*+ -> d+vt ) 

are known Q. Taking the R° = 0.616 and ii+ = 0.036 @, we obtain two 
sets of solutions for \X'/g\ and \X/g\, which gives two solutions for |A'/A|. The 
first is |A'/A| = 0.77 and the second is |A'/A| = 0.21 @. Taking A = -0.34 we 
get four possibilities for A' = ±0.26, ±0.071, which all have to be considered. 

III. THE AMPLITUDES AND BRANCHING RATIOS 
FOR NINE D -> Vl+l~ DECAYS 

A. The amplitudes 

In this section we turn to the amplitudes and branching ratios for the nine 
D -» Vl + r decays. The interaction Lagrangians (|), (H), (||), ©, (||) and 
the weak currents (0), (18) provide us with the vertices in the kinematical 



region, where the heavy quark and chiral symmetry are good (i.e. the velocity 
of the heavy mesons changes only slightly in the interaction and the energy of 
the light mesons is small). The problem is how to extrapolate the amplitudes 
to the rest of the kinematical region allowed in D — > Vl + l~ decays. We assume, 
that the vertices do not change significantly throughout the kinematical region, 
which is a reasonable assumption in D decays. At the same time we use the 
full heavy meson propagators l/(p 2 — m 2 ) instead of the HQET propagators 
l/(2mv ■ k). We account for the short life time of the intermediate neutral 
vector meson Vq by using the Breit Wigner form for the Vq propagator 

„• Yq 



q 2 - m\r o + iT Vo m Vo 

where Ty is the decay width of the Vo meson and q is its momentum. Then, 
using the interaction Lagrangians ([[]), ©, (|), (0), (0), @, (0) and the 
weak currents (|l^) , (|l8| ) , the calculation of the amplitudes for long distance di- 
agrams on Fig. 1 and short distance diagrams on Fig. 2 is straightforward. The 
calculated amplitudes for different diagrams have in general different Lorentz 
structure. It is convenient to treat the amplitudes of similar Lorentz structure 
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together. Then, the sum of the amplitudes within the model described above 
is given by the expression 

A [D{p) -> V(e {v) ,p {v) ) l + (p + ) r(p_)] = -%e 2 f Cab ~ 2 (19) 

V 2 q 

x e (y)/3 u(p_) 7l/ v(p + ) [e^ a/3 ^p Q ^p C + %A%\ , 

where q = p_ + p+ is the momentum of the intermediate virtual photon and 
the corresponding Cabibbo factors fc a b are given in Table 1. Ape and A p v v 
correspond to parity conserving and parity violating amplitudes, respectively. 
They get contributions from different diagrams in Fig. 1 (long distance) and 
Fig. 2 (short distance). The short distance amplitudes Ape and A p u v for 
the Cabibbo suppressed decays are given in Appendix A. The long distance 
amplitudes Ape and A p V y for the Cabibbo allowed, suppressed and doubly 
suppressed decays are given in Appendix B. 

The decay width for D — ► Vl + l~ is given by the square of the amplitude, 
summed over the polarizations of the three particles in the final state and 
integrated over the three body phase space 

r 1 \- [ \ M m2 d 3 p {v) d? P+ d?p- 

r = 2^(2^ £ J \^iv)^v-)\ l^^r^r SiPiv^+P-p) • 

(20) 



B. Discussion of the results 

Firstly, we present the results for the decays with the muon final state 
D — * V and we comment on the decays D — * Ve + e~ in the end. The 

branching ratios for the Cabibbo allowed, suppressed and doubly suppressed 
D — > Vn + fi~ decays are presented in Table 1. The last column presents the 
experimental upper bounds 14|. The other columns present our theoretical 



predictions, where the error bars are due to the uncertainty of the model pa- 
rameters A' and Cyvn, which can have any of the values A' = ±0.07, ±0.26 and 
Cvvn = ±0.31. The total branching ratio Br (total) containing long (Fig. 1) 
and short distance (Fig. 2) contributions is given in the fifth column. The third 
column presents the short distance part of the branching ratio Br(SD) calcu- 
lated from (H), which is present only in the Cabibbo suppressed decays. The 
fourth column presents only the nonresonant part of the long distance contribu- 
tion Br(LD 

nonr)- This part is bigger for the charged D meson decays, where it 
is mainly due to the diagram IV of Fig. 1. For the neutral D meson decays, the 
diagram IV vanishes and the remaining nonresonant diagram /// has smaller 
amplitude, which is proportional to A'. The parameter A' = ±0.07, ±0.26 has 
large uncertainty and we are only able to quote the upper limit for Br(LD nonr ). 
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Apart from the Cabibbo structure, the branching ratios depend mainly on 
whether the initial (and final) state is charged or neutral, with bigger branching 
ratio in the former case. We present also the distributions (1/Td)(1T(D — > 
V ii + / dq 2 as a function of q 2 (q 2 is invariant mass) for the typical 

representatives of the Cabibbo allowed (L>+ — ► in Fig. 3 and D° — > 

K*°[i + fi~ in Fig. 4) and suppressed (D° — > p fj, + fj,~ in Fig. 5 and Z?+ — > 
K* + p + p~ in Fig. 6), neutral or charged D meson decays. The short distance 
contribution (dot-dashed line) due to c —* ul + l~ transition is present only in the 
Cabibbo suppressed decays and it turns out to be much smaller than the long 
distance contribution. Concerning the long distance contribution, the resonant 
part is bigger than the nonresonant part (dashed line), except perhaps in the 
case of charged D meson decays at the low q 2 (see Figs. 3-6). Note, that 
the nonresonant LD contribution is generally smaller than the resonant LD 
contributions even in the regions well outside the resonance peak at q 2 = my - 
It is interesting to remark, that the short distance and nonresonant long distance 
contributions are comparable for Cabibbo suppressed neutral D meson decays 
D° -> p°l+r and D° -» ujI+1'. 

In the decays with the electron final state D — ► Ve + e~ the lowest kinemat- 
ically allowed q 2 is = (2m e ) 2 , which is smaller than g^ in = (2m At ) 2 in the 
D — ► Vfi + [i~ case. In the region q 2 > (2m fJi ) 2 the electron rates are practically 
equal to the muon rates. In the region q 2 < (2m^) 2 , however, the rates for 
D — ► Ve + e~ are extremely enhanced due to the photon propagator 1/q 2 . How- 
ever, the region down to <2^ n = (2m e ) 2 requires a more accurate treatment of 
the q 2 dependence when q 2 approaches to , which is beyond our scope here. 
We have calculated the D — * Ve + e~ branching ratios with the lower cut off 
q 2 = (2m^) 2 and have obtained values which are very close to the D — ► Vfi + [i~ 
branching ratios (the D — ► Vfi + fi~ branching ratios are obtained integrating 
over the whole q 2 = [(2m At ) 2 , (mo — my) 2 ] region ). 

The Cabibbo allowed decays D° — > K*°/i + fj 1 ~ and — > p + fi + fi~ with the 
predicted branching ratios of the order 10 -6 and 10 -5 , respectively, have the 
best probability for their early detection. Note that their branching ratios are 
not far below the present experimental upper bound. 

In the Cabibbo suppressed decays, the short distance contribution due to 
FCNC transition c — > ul + l~ has branching ratio of order Br(SD) ~ 10 -10 and 
is therefore well masked by the long distance branching ratios of order 10 -7 . 
Obviously, to observe the FCNC transition c — > u within the Standard Model, 
one must most likely look for other possibilities. Still, new physics could en- 
hance the SD part to be of the same order as the LD part or bigger |^, [|, |j . In 
this case the branching ratios well above 10~ 7 for Cabibbo suppressed decays 
D — > V/j, + fi~ would signal new physics. As the present experimental upper 
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bound is much higher, these decays still contain a large discovery window. 
IV. SUMMARY 

We have calculated the long and short distance contributions for nine D — > 
Vl + l~ decays within the Standard Model. The short distance contribution is 
present only in the Cabibbo suppressed decays and is due to the flavour chang- 
ing neutral transition c — ► ul + l~ . The long distance contribution is composed 
of the resonant part, which arises from the intermediate light vector meson Vq 
exchange (D — ► VVq — > V*y — > V7 + Z _ ), and the nonresonant part, which arises 
from the direct photon emission {D — > — > The branching ratios are 

calculated using an effective Lagrangian, which combines heavy quark symme- 
try and chiral perturbation theory, and are given in Table 1. The most frequent 
decays are the Cabibbo allowed decays, which are expected at the rates, that are 
not much lower than the present experimental upper limit: Df — ► p + p + p~ is 
expected at the branching ratio of approximately 3-10 -5 , while D° — > K*°fj, + fj,~ 
is expected at 1.7 • 10~ 6 . The Cabibbo suppressed decays on the other hand, 
are typically expected at [3 — 7] • 10~ 7 range for D° — ► p° (uj) fi + /j," and Df — > 
K* + p + p~ decays and in the 10~ 6 range for D + — > p + p + p~ decay. Accord- 
ingly, branching ratios well above 10 -7 for Cabibbo suppressed decays could 
signal new physics. In all the Cabibbo suppressed decays the short distance 
contribution is well masked in the Standard Model by the resonant long dis- 
tance contribution. In the case of D° — > p°(u)l + l~ decays, however, the short 
distance contribution is of comparable size as the nonresonant long distance 
part. 
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fCab 
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< 1.9 icr 8 


[i fi - 1 91 10~ 6 
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a\(? 





4.0 10~ 6 


[3.0 - 3.3] 10~ 5 




L> u -> p u // + /* _ 


—a<isc 


9.7 l(T iU 


< 4.8 l(T iU 


[3.5 - 4.7] 10~ y 


< 2.3 10~ 4 




— Q^SC 


9.1 10~ 1U 


< 3.7 10~ 1U 


[3.3 - 4.5] 10- y 


< 8.3 10~ 4 









< 1.1 1(T 9 


[6.5 - 9.0] 10~ 8 


< 4.1 10~ 4 


D+ -> p + p, + p~ 


— a\sc 


4.8 l(T y 


2.7 10"' 


[1.5 - 1.8] 10 -e 


< 5.6 10~ 4 


D + _ K *+p+p,~ 


a\sc 


1.6 10~ y 


1.5 10-' 


[5.0 - 7.0] 10~ 7 


< 1.4 10 :! 


D+ -» K*+p+p~ 


— ais 





1.0 10~ 8 


[3.1 - 3.7] 10~ 8 


< 8.5 10~ 4 


L> u -» K^p+p' 


-a 2 s 2 





< 5.0 lO" 11 


[4.4-5.1] 10~ 9 





Table 1: The branching ratios for the Cabibbo allowed, suppressed and doubly 
suppressed D — > Vp + p~ decays. The last column presents the experimental 
upper bounds |l3|, p^] , while the other columns present our theoretical predic- 
tions. The total branching ratio Br{total) containing long (Fig. 1) and short 
distance (Fig. 2) contributions is given in the fifth column. The third col- 
umn presents only the short distance part of the branching ratio Br(SD). The 
fourth column presents only the nonresonant part of the long distance contri- 
bution Br(LD nonr ). The error bars in the Table are due to the uncertainty 
of the model parameters expressed by the possibilities A' = ±0.07, ±0.26 and 
Cvvu = ±0.31. The branching ratios for D — > Ve + e~ obtained with the lower 
cut off q 2 = (2m M ) 2 (q 2 the invariant p + p~ mass) are almost exactly the same as 
the branching ratios for D — ► V p + p~ given in this Table. The second column 
gives the corresponding Cabibbo factors fc a b in terms of the Cabibbo angle 
c = cos 0c and s = sin^c". 



APPENDIX A: The short distance amplitudes 



In this Appendix we list the values of the Willson coefficients Aj, Bi and F\ 
in the short distance Lagrangian Csd (§) and give the resulting short distance 
amplitudes Ape and A?p V from (|lS|). 

The coefficients A> Bi and F\ have been obtained in the leading order by 
Inami and Lim 2J] and following the notation of @ one has 



A 

B; 



(21) 



c box + C Z _ sin 2 0w ( F i + C Z} 

- sin 2 OwiFt + Cf) , 

where C\ ox Cf , F\ and F\ are kinematic factors, which depend on the ith- 
quark mass through = mf j 

x j In x- j 



a 



box 



1 



Cf 



T 



+ 



8 Xi 



+ 



1) 

3 2x 2 



22 



X i 



(xi - If 



Inxi + j(£,Xi 



13 



+ -x 



+ X 




3 xf In Xj 



-l) 2 ' 2 (Z 

The summation i in Csd (||) runs over down-like quarks (d, s and b) to which 
charm can couple, while Q = — 1/3 is the corresponding charge of the inter- 
mediate quarks (we note that F\ and F2 have been calculated in Q using the 
wrong charge Q = 2/3). The gauge dependent term 7^, xi) p8[ cancels out in 



(23) 



the combinations A4, B{ and F\ (21). Since the ratios xj, x s and Xf, are of orders 
1CP 8 , 10 -6 and 10~ 3 respectively, the terms proportional to the powers of x% 
can be safely neglected in (22). With this approximation C\ ox = Cf = —3/8, 
F[ = -2 In 2*/ (9a* 
occurs and we obtain 



9) and F% vanishes. In this limit the GIM cancelation 



Y,V i A i = Y,V i B i = A 



SD 



-0.065 



(24) 



Consequently, the short distance Lagrangian (j||) effectively contains only the 
vector lepton current Ij^l but not the axial vector I'j^'jsl one. 

The short distance c — > ul + l~ transition contributes only to the Cabibbo 
suppressed decays. Here we give the short distance contributions for the parity 
conserving Ape and parity violating A p u v amplitudes, which are needed to cal- 
culate the Cabibbo suppressed amplitudes A[D(p) — > V(e(y),P(y)) l + {p + ) Z - (p_)], 
(|l9|). Within the model used, these amplitudes are given by the diagrams on 
Fig. 2: 



A P c(D ^p°l+r) 
A PC (Df^K*+l + r)/V2 
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A p v c {D Q ^cjPin-) = A%(D 



, 9V ym D 
fc'ab 

Pv I n0 



ui + r) = 



a 2 - 



qV 1 



m 



D 



-r) = . 



-^p + l + l-)/V2 

A SD g 2 
16vr 2 sin 2 8 W ' 

(25) 
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The last equation is a result of Eq. (||) and the quark content of the <j> meson. 
The relevant constants are presented in Appendix B. 



APPENDIX B: The long distance amplitudes 



In this Appendix we give the expressions for the parity conserving Ape and 
parity violating amplitudes Apy, which are needed to calculate the amplitudes 
A[D(p) -> V(e {v) ,p (v) ) l + (p+) r(p_)] dH) for nine D -» Vl + l~ decays. The 
following amplitudes Ape and Apy contain the long distance resonant and 
nonresonant contributions coming from the diagrams on Fig. 1. The coefficients 
and constants needed for the evaluation of the amplitudes will be given bellow: 
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(26) 



Here q = p~ + p + and my can be approximately taken as the average of the 
6, uj and p masses. The coefficients J D , K v , L v , M D and N are expressed as 
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and 

N = 9p gl H _ 

q 2 - m 2 + iT p m p 3(q 2 - m 2 + i^m^) 3(<? 2 -m^ + iT^m^) 

The functions /i mix , f' lmix , hmix and are denned by 

/„ 1 + c 2 sc _ sc 1 + s 2 

r l - 5c 2 5sc , __ ,/y -5sc 1 - 5s 2 

V8 Jn J v > V8 Jr? jy 

where s = sin#p, c = cos9p and #p ~ 20° is the r\ — rj mixing angle. The 
values of the masses, decay constants and decay widths used are given in Table 
2. 
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4> 


1.02 
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Table 2: The pole masses m, decay constants / and decay widths T in GeV; 
the constants gy in GeV 2 . 
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Figure Captions 

Fig. 1. Skeleton diagrams of various long distance contributions to the de- 
cay D — > resulting from Eq. The spectator diagrams of type Aspec,y 
(see Eq. (§)) are shown in Fig. la, the spectator diagrams of type As pe c,v 
are shown in Fig. lb and the weak annihilation diag rams are shown in 
Fig. lc. Different diagrams are denoted by the roman numbers I — VIII. The 
square in each diagram denotes the weak transition due to the long distance 
Lagrangian Cld (Eq. @). This La grangian contains a product of two weak 
currents, each denoted by a dot. 

Fig. 2. Skeleton diagrams of short distance contributions to the decay 
D — > Vl + l~ due to c — > ul + l~ transition. The squares in the diagrams denote 
the weak transition due to the short distance Lagrangian Csd (Eq. |||). This 
Lagrangian contains a product of a quark and a lepton weak currents, each 
denoted by a dot. 

Fig. 3. The differential branching ratio (l/Y D )dT(D° -> K*°n + p-)/dq 2 
as a function of q 2 (q 2 the invariant ft + /i~ mass). The full line corresponds 
to the total branching ratio, while the dashed line represents the nonresonant 
long distance part. In the calculation the model parameters A' = 0.26 and 
Cvvu = 0.31 were used. 

Fig. 4. The differential branching ratio (l/Ti))dT(Df — ► p + p + /i~)/dq 2 
as a function of q 2 (q 2 the invariant mass). The full line corresponds 

to the total branching ratio, while the dashed line represents the nonresonant 
long distance part. In the calculation the model parameters A' = 0.26 and 
Cvvu = 0.31 were used. 

Fig. 5 The differential branching ratio (l/T]j)dT(D — > p°p + p~)/dq 2 as 
a function of q 2 (q 2 the invariant /J> + i-i~ mass). The full line represents to the 
total branching ratio, the dot-dashed line represents the short distance part, 
while the dashed line represents the nonresonant long distance part. In the 
calculation the model parameters A' = 0.26 and Cwn = 0.31 were used. 

Fig. 6 The differential branching ratio (l/T D )dT(D+ -> K*+ p+ yr)jdq 2 
as a function of q 2 (q 2 the invariant mass). The full line represents to the 

total branching ratio, the dot-dashed line represents the short distance part, 
while the dashed line represents the nonresonant long distance part. In the 
calculation the model parameters A' = 0.26 and Cvvu = 0.31 were used. 
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